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The choice of the 4-velocity in the relativistic hydrodynamics proposed by Landau and Lifshitz 
is demonstrated to be the only one satisfying all general principles. Especially the far more widely 
adopted Eckart choice has to be abandoned. 
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Three reasons come to mind for seeking a covariant 
formulation of the hydrodynamic theory for simple flu- 
ids: In a fundamental vein, we want hydrodynamics as a 
basic theory to be covariant; in fact, this was one of the 
first few tasks tackled after the birth of special relativity. 
More practically, relativistic hydrodynamics is increas- 
ingly employed in cosmology and astrophysics to study 
dissipative processes, such as the relaxation of inhomo- 
geneities in the early universe, or its viscosity-driven, in- 
flationary expansion |jj . Finally, any residual ambiguities 
in the relativistic hydrodynamics are an indication of our 
less-than-perfect grasp also of the Galilean version. They 
require close scrutiny. 

Covariant hydrodynamics can be found in most text- 
books on relativity and astrophysics. Better ones give 
two versions, one due to Eckart ||||, the other by Lan- 
dau and Lifshitz (LL) While Eckart pins the macro- 
scopic 4-velocity to the total particle current, LL set it 
proportional to the total energy flux. So there is no dis- 
sipative particle current in the Eckart version, in direct 
analogy to the non-relativistic case; while LL — at first 
sight somewhat odd — do have one, and they refrain from 
a dissipative energy current instead. In the literature, the 
Eckart version is much more widely employed 0J^,|^1, it 
seems the more traditional theory. Yet, both reduce to 
the familiar non-relativistic hydrodynamics for c — » oo; 
besides, standard textbooks (eg. Weinberg ||) regard 
the two versions as equivalent, as being related by a sim- 
ple transformation of the velocity. Clearly, one need not 
worry about a mere difference in the reference frame. 

This equivalence, however, is a fallacy. Let us recall 
how the velocity is defined in hydrodynamic theories, 
consider first the non-relativistic case. The standard 
Gibbs relation, valid not only when equilibrium reigns, 
takes the entropy density s as a function of five conserved 
densities, energy e, momentum g, and mass g, 

T ds = de — ~v-dg — [idg . (1) 

It states unequivocally that the velocity is a thermody- 
namic quantity, v = —T(ds/dg), known if the local en- 
tropy density is. It contains only equilibrium informa- 
tion. This is of course the concept of local equilibrium, 
one of the few founding principles of the hydrodynamic 



theory: It takes far less time to establish equilibrium lo- 
cally than globally; the first is a microscopic time r (re- 
ferred to as the collision time in dilute systems) , the other 
grows with the system dimension and is macroscopic. As 
long as the frequency is small enough, wt < 1, the Gibbs 
relation holds, and all thermodynamic variables e, g, g, 
and their conjugate variables T, v, /i contain only in- 
formation about the local equilibrium state. Especially, 
they possess a well defined parity under time reversal. 
The relativistic description is hardly different: The ve- 
locity v becomes the 4-velocity, = —T(ds v /dT^) : 
with and T^ v being the equilibrium entropy 4-vector 
and energy-momentum tensor, respectively. So too 
contains only equilibrium information. 

Now, the difference between the two 4- velocities men- 
tioned above is clearly in dissipative quantities, so at 
least one contains non-equilibrium information and can- 
not be the correct hydrodynamic velocity. This argument 
agrees with the prima facie evidence that the two ver- 
sions of relativistic hydrodynamics have different types 
of differential equations. For example, the equation for 
the transverse velocity is elliptic in the Eckart version, 
and parabolic in that of LL. (In the so-called extended 
thermodynamic theories — cf. discussion below — the 
equation can be rendered hyperbolic for either choice of 
the 4-velocity. Still, as these are extensions of different 
hydrodynamic theories, they remain distinctly different.) 

Having clarified that both versions are inequivalent, we 
obviously need to address the question: Which, if any, is 
the correct theory for relativistic dissipative fluids? In 
seeking uniqueness in a relativistic theory — since the 
c — > oo limit proves inconclusive — it is natural to ex- 
amine how uniqueness is achieved in the Galilean ver- 
sion. Unfortunately, this is something of a red herring, 
as the lack of a dissipative particle current in the Galilean 
hydrodynamics, or j = qv, is more a statement of mi- 
croscopic plausibility, and maybe the bold summary of 
countless experiments; it is not the result of a cogent and 
general deduction. In fact, a classic paper by Dzyaloshin- 
skii and Volovik [|| proposes to include the dissipative 
term j — g v ~ V/i. 

Nevertheless, a footnote by LL, in §49 of [||, states 
that the mass current j should always be equal to the 



1 



momentum density g = gv, and therefore cannot pos- 
sess any dissipative terms. Their line of argument, re- 
grettably, falls short of being ironclad: Starting from the 
continuity equation and the center-of-mass motion, they 
assert the validity of J dV g = J dV ], where the integra- 
tion volume is that of the total system. The reader is left 
wondering especially about the alleged equivalence to the 
local relation g = j. 

In this paper, we shall ameliorate both aspects. We 
provide a clear cut proof of g = j, demonstrating the 
rigorous validity of the standard form of the Galilean 
hydrodynamics; and we demonstrate that only the LL 
version of the relativistic hydrodynamic theory conforms 
to all general principles. The proof g — j takes place in 
very much the same way as deducing the symmetry of 
the momentum flux, or stress tensor, from local angular 
momentum conservation. The relevant conserved den- 
sity here is i = g x — g t, with x and t denoting the space 
and time coordinates, respectively. While i is known to 
be additive, conserved, and the direct consequence of the 
invariance under Galilean boosts , it has not been hith- 
erto included in thermodynamic and hydrodynamic con- 
siderations. As will become clear soon, this is a serious 
omission. Its inclusion not only establishes the form of 
the mass flux, but also leads to general thermodynamic 
equilibrium conditions that are valid for any reference 
frame. Surprisingly, these simple yet fairly fundamental 
relations are new. 

Relativistically, the information of g = j is automati- 
cally included in the symmetry of the energy-momentum 
tensor, though j is now the inertial mass current, ie. the 
total energy flux including the rest mass. Since the mo- 
mentum density remains a thermodynamic variable, with 
a negative parity under time inversion, neither the mo- 
mentum density nor the inertial mass current may con- 
tain any dissipative terms (in the local rest-frame) . This 
excludes any covariant theory that does not adopt the 
LL 4- velocity. 

Two pieces of information were needed in each of the 
above cases: (i) The equality of the momentum den- 
sity to the (rest or inertial) mass current; (ii) the fact 
that the momentum density is a thermodynamic vari- 
able with a well defined parity. It is ironic that while the 
condensed matter people were wondering about the first 
piece, which for the relativists is a trivial consequence of 
the 4-notation, the latter ignored the second, something 
the former group never does. We physicists are indeed 
a community divided by notations. Piercing both pieces 
together, the kinship to be expected between the two ver- 
sions of hydrodynamics becomes evident. For instance, 
repeating the relativistic mistake in the non-relativistic 
theory, ie. violating the second condition while upholding 
the first, leads to, as will become clear soon, 

g=i = gv-x(Vfi + d t v). (2) 

This is in striking similarity to the momentum density in 



the Eckart theory, the dissipative part of which has the 
form r l4 = -x{T- x diT + d t Ui) f§. Yet Eq. §) is man- 
ifestly unphysical: The total and conserved momentum 
must remain J d 3 x gv, irrespective whether the system is 
in equilibrium or not, or what its acceleration is. 

Of the three issues plaguing the relativistic hydrody- 
namics — uniqueness, causality and stability — we focus 
on the first. But we need to comment on the other two, 
as they have been the starting points of worthwhile ef- 
forts in the past that partially tic in with our results. 
First, causality. Strictly speaking, the diffusion equation 
implies signals with infinite velocity, or horizontal "world 
lines". While unphysical generally, this defect is aggra- 
vated in relativity: When viewed from a different frame, 
the world lines tilt, implying signals that go backwards 
in time. To repair this, extended thermodynamic theo- 
ries [^]]|] were put forward which start from the hydrody- 
namic theories but include additional dynamic variables. 
The resultant larger set of coefficients can be chosen such 
that all the differential equations are hyperbolic, ensur- 
ing causality. The price for this nice feature is a rather 
more complicated theory, and the difficulty of finding a 
universally valid and accepted set of additional variables 
— except perhaps in dilute systems. 

But we may also take a more perspective view, and 
accept that the diffusion equation is not an exact math- 
ematical statement. Rather, it is an approximative de- 
scription — confined to the hydrodynamic regime, with 
an accuracy of the order of thermal fluctuations. Tak- 
ing this into account, (eg. considering only amplitudes 
of the variables that are above a minimal threshold,) the 
signal velocity never exceeds that of the constituent par- 
ticles Pj, excluding any acausal consequences. 

Next, stability, first in the fluid's rest frame: The LL 
theory is stable with respect to small fluctuations around 
an equilibrium configuration, not so the Eckart version 
fiof , and remarkably, nor the non-relativistic theory that 
contains Eq. (||). In fact, both suffer from the same prob- 
lem. Consider a small but spatially homogeneous veloc- 
ity field with Vp, V/i = 0, the Navier-Stokes equation 
reduces to dtg = 0, or gdtv — x<9 t 2 v = 0, which (in ad- 
dition to the usual v = const) obviously also contains 
the run-away solution ~ e^slx)*- _ Similarly, with a mo- 
mentum density that contains the acceleration dtUi, the 
Eckart choice cannot help to avoid an analogous instable 
solution. 

This would represent an independent argument favor- 
ing the LL choice, except that — as observed by Hiscock 
and Lindblom — in frames moving with respect to 
the fluid the diffusion equations in the LL theory also de- 
velop diverging solutions, which grow exponentially with 
microscopically short characteristic times. For lack of 
space, we briefly summerize our reasons for believing that 
this frame-dependent instability does not constitute suf- 
ficient ground to reject the LL choice, and promise a 
detailed account in a forthcoming paper. Consider the 
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parabolic diffusion equation, dt'd — ad^d = 0. Its char- 
acteristics are the lines t = const, and only if the initial 
values are prescribed on one of these, do we have a sin- 
gle bound mode, 6fl e ikx ~ ak2 \ with e K @. Initial 
data on a non-characteristic curve, say x + (3 t — const, 
(3 e H, generally produce two independent solutions. For 
the simplest case of /3 = 0, they are S'&i e -' l ( nx +^ t ) e nx 
for x < and 8d 2 e^-^ e~ nx for x > 0, with = 
(w/2a) 1 / 2 , u> £ TR. In the respective wrong region, one 
solution appears unbound. Being invariant with respect 
to coordinate transformations, the characteristics of the 
boosted diffusion equation that Hiscock and Lindblom 
consider are t — j(t Isi) = const, with t the proper 
time. The solutions they examine, however, satisfy ini- 
tial data on the non-characteristic t = const, where t 
is the time in the moving frame. So the appearance of 
an unbound solution for t — > oo is a mathematical conse- 
quence to be expected. Nevertheless, the diverging mode, 
being absent in the rest frame, must not be observable 
in a moving one. And it is not, as it only exists for 
negative times t, and decays for t — > — oo within a micro- 
scopically brief period that is outside the hydrodynamic 
regime. In fact, this mode is just one of those signals 
discussed above that run backwards in time in moving 
frames. (These arguments do not apply to the Eckart 
instability. It happens in the fluid's rest frame, where 
any deviation from the non-relativistic hydrodynamics is 
worrisome.) 

The extended theories are stable for both choices of the 
4-velocity if linearized; though the Eckart version turns 
instable again if non- linear terms are included p"2| . 

We conclude: Within its range of validity, the rel- 
ativistic hydrodynamics is just as healthy as the non- 
relativistic theory. If someone is willing to put up with a 
few acausal consequences, blatant but recognizably out- 
side this range, he retains the benefit of a simpler theory. 
If not, he may resort to the extended theory — though 
it has to be one that adheres to the LL choice of the 
4-velocity. 

Let us now consider the hydrodynamics in greater de- 
tails, starting again with the non-relativistic version. The 
equations of motion for the thermodynamic variables of 
Eq. © are, 

d t s + Vf = R/T, 9ie + V-q = 0, (3) 

d t9l + d k n lk = o , d tQ + v-j = o. (4) 

We explicitly include the conserved quantity 

1 = fd 3 x(gx-gt) = MX(t)-Gt (5) 

in our consideration, where M , G, and X(t) denote the 
total mass, the total momentum, and the center-of-mass 
coordinate, respectively. Clearly, I/M is the initial coor- 
dinate of the center of mass, so we may perhaps refer to 
I as the center-of-mass inertial coordinate (comic), and 
to i = g x — g t as the COMIC density. 



Neither the angular momentum nor the COMIC re- 
quires an independent equation of motion. Writing 

^t{^ikm %k 9m) — dni^ikm limn) ~t~ ^ikm^-mk: One 

finds that the angular momentum density £ikm %k 9m 
obeys a continuity equation only if IL,fc = Hki- Analo- 
gously, d t (gxi -gtt) = -d k (j k %i - Ihk t) + ji - g { holds 
for the COMIC density, a locally conserved quantity, hence 
g = j. This concludes the clear cut and simple proof we 
were looking for. 

Next we deduce thermodynamic equilibrium condi- 
tions including the conservation of the COMIC I. Max- 
imizing the total entropy S = J d 3 x s subject to the con- 
servation of energy, momentum, mass, angular momen- 
tum, and COMIC, we have J d 3 x {5s — Ai 5e + Ai • Sg + 
A2 Sg + A 2 -S(x x g) — A 3 -S(gx — gt)} — 0, where the 
eleven coefficients Ai^ and Ai 2,3 are constant Lagrange 
parameters. Employing Eq. (|lj), we deduce, for arbitrary 
variations 5e, 5g and Sg (with Sx = St = 0) , 

1/T = X 1 , n/T = A 2 - A 3 -x, (6a) 
v/T = A x + A 2 x x + A 3 t. (6b) 

The last expression does not imply an accelerating mo- 
mentum, as Dixon concluded in Ch.4 §4d of Ref. ||. To 
see this directly, consider uniform space and time trans- 
lations: Setting now Sx, St = const, and requiring that 
the equilibrium conditions remain unaltered, we arrive at 

M A 3 = -A 2 x G (6c) 

and the dependent A3-G = 0. Together, the equi- 
librium conditions (Q) are explicitly Galilean covariant: 
Introducing the chemical potential fio = u + \v 2 of 
the local rest-frame, they can be expressed as T = T, 
Mo = M + x (x - X)] 2 , v = V + Q x (x - X), 
with X = X(t) being the center-of-mass coordinate, and 
T, Ji, V, $1 redefined constants. Clearly, /j, (and hence 
the density g) only depends on the rotation velocity in 
the center-of-mass frame, and not on the center-of-mass 
motion. Without including the COMIC I (ie. setting 
A3 = above) LL obtained, in a similar calculation |fl3| , 
v = V + fl x x, and concluded that the equilibrium ve- 
locity v of a general frame has to be a constant of time. 
But this is clearly only correct in special frames, when 

v||n. 

Now Eq. (||) is derived. First we remark that the un- 
usual form of the thermodynamic force V/i + <9 t v is a 
natural consequence of Eqs. (^|): This combination van- 
ishes in equilibrium and may therefore serve as a le- 
gitimate thermodynamic force. More technically, given 
the existence of j D = j — gv, both g and e acquire a 
dissipative part, g = g Eq + j D and e = e Eq + v • j D . 
Substituting Eqs. (^) and ^ into the Gibbs relation, 
T dts = dt£ Eq — v ■ dtg Eq ~ ndtg, one obtains the entropy 
production R = —j D - (V/i+9tv) + from which Eq. (Q) 
[with x > 0, such that R > 0] results. 
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Turning our attention to special relativity, the relevant 
hydrodynamic equations ([j]), (||), and (0) generalize to 

7W 1 = -u v dT v » - fidn^ , (7) 
d ll (s»+a») = R/T, (8) 
d u (T^ + t^) = , c^K + v") = . (9) 

Notations: Greek indices run from 1 to 4, Latin indices 
only to 3; the speed of light is unity; the metric is rf v = 
diag(l, 1, 1, —1); the coordinate 4-vector is x M = (x, t), 
so dfj, = (V, dt)\ the 4-velocity is — j(v, 1) with 
7 = (1 — u 2 ) -1 / 2 , hence u M u M = —1; ct m , t^, and v* 1 
are the respective dissipative parts of the entropy 4-flux, 
energy-momentum tensor, and particle 4-flux, they have 
a different parity under time reversal from their reactive 
counterparts and vanish in equilibrium. 

In the laboratory frame, in which the local fluid veloc- 
ity is v, the reactive, equilibrium terms are: = su^, 
T^ v = (e + p)u^u v + prf- v ' , and = nu^, where 
p = —e + Ts + fin is the pressure, e = eo + Q the den- 
sity of total energy, n = g/m the particle number den- 
sity, and s as before the entropy density, all taken from 
a local comoving frame. The last three are related by 
Tds = de — fidn, so \i = m(deo/dg) + m has a differ- 
ent definition than in the Galilean case, but p does not. 
The relativistic version of the Gibbs relation, Eq. (]?]), 
is obtained by multiplying the rest-frame relation with 
it M . The conservation of the 4-angular momentum is en- 
sured by the symmetry of the energy-momentum tensor: 
Tt iu _ TVIi and rM „ = T „ M _ Ag discussedi this i nc l u des 

the equality of the momentum density and the total en- 
ergy flux. 

Now consider the explicit form of the dissipative terms 
a^, , and z/'L They are determined by the rate of en- 
tropy production R, a positive Lorentz scalar. Inserting 
Eqs. (^) and (||) in ([?]), and requiring R to be a sum of 
products of thermodynamic fluxes and forces, we arrive 
at TV 1 = -u v - n and 

R/T = -t^ d [p [u v) /T] - ^d^/T) , (10) 

where the bracket denotes symmetrization, eg. d^u^j = 
(d^Uv + d v Ufj)/2. Global equilibrium conditions are met 
if the two forces di^ [u v \ /T] and d^(n/T) vanish. Ir- 
reversible thermodynamics generally prescribes the On- 
sager ansatz, setting t^ v and as linear combinations of 
d(fj\u v )/T] and d^n/T) — subject to the requirements 
that thermodynamic variables do not possess any dissi- 
pative counterparts in the local rest-frame (lrf). The 
thermodynamic variables of Eq. (j?]) reduce to s 4 = s, 
= (g, e), and n 4 = n in the lrf. So the lack of 
dissipative counterparts implies cr 4 = t^ 4 = v A = 0, of 
which the covariant expressions are, 

UpO* =U V T»" = =0. (11) 

These are the conditions implemented by LL, and the 
ones we need to heed while evaluating Eq. (|Tc|). 



It must have been a source of confusion that g itself 
vanishes in the lrf — reducing T m4 to T 44 = e, and 
seemingly leaving only r 44 = 0. This overinterprets the 
lrf. What we actually need is to examine the infinites- 
imal changes of the variables, dT^ 4 = (dg, de), and un- 
derstand that the lrf does not imply dg = 0, as we must 
allow for dtg, digk ^ 0. Non-relativistically, of course, g 
being a thermodynamic variable is never disputed. 

It is not incidental that the conditions (^lj) rule out 
any time derivative in R. Violating (^lj), we find from 
Eq. ©: t* 4 ~ d t u u r 44 ~ d t T, v 4 ~ d t (p/T) in the 
lrf. In the equations of motion, each yields its own 
run-away solution, altogether five. Eckart's conditions, 
it M cr M = u M u v t^ v = = 0, partly violate Eqs. (|Ti"l), so 
it is not surprising that his "momentum of heat" r l4 ~ 
dtUi gives rise to the instable solution discussed above. 

The covariant hydrodynamic theory that entails the 
4-velocity as proposed by LL is the appropriate theory 
to employ if the velocity difference in the system is no 
longer small when compared to light velocity. However, 
it does not consider charges and electric currents that are 
frequently present in astrophysical systems. For this, one 
needs in addition the covariant version of the dissipative 
Maxwell equations H], to be published elsewhere. 
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